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INTERLAMINAR  STRESS  SINGULARITIES  AT  A  STRAIGHT 
FREE  EDGE  IN  COMPOSITE  LAMINATES 

* 

I.  S.  Raju  and  John  H.  Crews,  Jr. 

NASA  Langley  Research  Center 
Hampton,  Virginia  23665 

SUMMARY 

A  quasi-three-dimensional  finite-element  analysis  was  used  to  analyze 
the  edge-stress  problem  in  four-ply,  composite  laminates.  The  seven  laminates 
that  were  considered  belong  to  the  laminate  family  [6/ (0-90) ] s ,  where 
0  <  6  <  90.  Systematic  convergence  studies  were  made  to  explore  the  existence 
of  stress  singularities  near  the  free  edge.  The  present  analysis  appears  to 
confirm  the  existence  of  stress  singularities  at  the  intersection  of  the  inter 
face  and  the  free  edge.  The  power  of  the  stress  singularity  was  the  same  for 
all  seven  laminates  considered. 


INTRODUCTION 

Many  composite  laminates  develop  high  interlaminar  stresses  near  free 
edges  due  to  mechanical  loads,  thermal  loads,  or  hygroscopic  effects.  The 
interlaminar  stresses  can  lead  to  delamination  and  eventual  failure  of  the 
laminate.  Several  numerical  analyses  have  computed  stress  distributions  in  a 
finite-width  laminate  subjected  to  uniform  axial  strain.  This  case  is  often 
referred  to  as  the  straight-edge  or  the  edge-stress  problem.  In  reference  1,  a 
finite-difference  scheme  was  used  to  analyze  the  edge-stress  problem  of  a  [±45] g 
laminate.  In  reference  2,  a  three-dimensional  finite-element  method  was  used  to 
analyze  the  edge  effects  in  [±45]s  and  [90/0] g  laminates.  In  reference  3,  a 
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quasi- three -dimensional  finite-element  analysis  was  used  to  analyze  [0/90]  , 

[90/0]  ,  [±45]  ,  and  several  quasi-isotropic  laminates.  These  analyses  used 

s  s 

coarse  mesh  models  because  of  limitations  on  computing  storages  and  times. 
However,  all  showed  that  the  interlaminar  stresses  had  very  localized  peaks 
near  the  free  edge.  If  these  localized  peaks  were  caused  by  a  singularity, 
as  suggested  in  references  1  and  3,  singularity  analyses  could  be"developed 
and  used  to  study  delaminations. 

The  present  study  investigated  the  existence  of  stress  singularities  at  a 
free  edge  in  several  four-ply  laminates  that  belong  to  the  laminate  family 
16/(6-90)]  ,  0  <  0  <  90.  To  this  end,  quasi-three-dimensional  finite-element 

analysis  was  used.  The  analysis,  which  is  similar  to  that  of  reference  3, 
used  eight-noded  isoparametric  elements. 

The  interlaminar  stresses  are  presented  as  distributions  through  the 
thickness  at  the  free  edge,  and  along  the  interface  between  the  plies.  The 
results  are  compared  and  discussed  for  seven  different  laminates.  Stress 
distributions  calculated  from  the  analyses  were  used  in  a  procedure  to 
investigate  the  existence  of  stress  singularities. 

SYMBOLS 

strength  of  singularity,  MPa 

half-width  of  the  straight-edge  laminate,  m 
Young’s  modulus  in  i-direction,  MPa 
shear  modulus,  MPa 
ply  thickness,  m 

nondimensional  distance  from  the  free  edge,  r  =  (b  **.  y)/h 
displacement  functions,  m 


E.  . 

li 

G.  . 
ij 

h 

r 

u,v,w 
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u,v,w  displacements  in  x-,  y-,  and  z-directions,  respectively,  m 

x,y,z  Cartesian  coordinates,  m 

a  power  of  the  singularity 

e  uniform  axial  strain  imposed  on  the  straight-edge  laminate, 

°  e  =  0,001 

o 

0  outer  ply  angle,  deg 

V.  Poisson's  ratios 

ij 

{a}  Cartesian  stresses,  O x#  a  ,  0^9  Cyz>  MPa 

Subscripts: 

i,j  iJ  =  1,2,3 

1,2,3  longitudinal,  transverse,  and  thickness  directions  of  a  unidirec¬ 

tional  ply 

ANALYSIS 

This  section  describes  the  laminates  analyzed,  the  quasi— three— dimensional 
formulation,  and  the  finite— element  idealizations  used  in  this  study.  This 
section  also  discusses  the  procedure  used  to  investigate  the  existence  of  a 
stress  singularity. 


Laminates 

The  laminates  considered  in  the  analyses  were  [0/90] g,  [15/-75]g, 

[30/-60]  ,  [ ±45]  ,  [60/-30]  ,  [75/-15]  ,  and  [90/0]  .  As  previously  mentioned, 

s  s  s  s  s 

these  seven  laminates  belong  to  the  family  [0/(0-9O)]g,  0  <  0  <  90.  The  top 

ply  angle,  -  0,  is  shown  in  figure  1(a). 

Each  ply  was  idealized  as  a  homogeneous,  elastic  orthotropic  material 

with  properties  (refs.  1,  2,  and  3), 
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=  137.9  GPa  (20  x  10  psi) 


E33  =  14.48  GPa  (2.1  x  10  psi) 


G23  =  G3  =  5.86  GPa  (0.85  x  10 


V13  "  V23  =  °'21 


psi) 


The  subscripts  1,  2,  and  3  correspond  to  the  longitudinal,  transverse,  and 
thickness  directions,  respectively,  of  the  individual  ply. 


Quasi-Three-Dimensional  Analysis 

Figure  1(a)  shows  a  long,  symmetric  laminate  loaded  in  the  x-dir ect ion . 
The  laminate  has  a  width  of  2b  and  has  four  plies,  each  of  thickness  h. 

Away  from  the  ends,  the  displacements  in  any  x  =  constant  plane  were  assumed 
to  be 

u(x,y,z)  =  eQX  +  U(y,z)  ' 

v(x,y,z)  =  V (y , z)  >  (1) 

w(x,y,z)  =  W(y, z)  J 

where  £  is  a  uniform  axial  strain,  and  U,V,W  are  functions  of  the  coor- 
o 

dinates  y  and  z  alone  (ref.  1).  Furthermore,  because  of  symmetries 
inherent  in  the  layups  considered,  U,  V,  and  W  satisfy  the  following 
requirements : 

U(y,z)  =  -U(-y,-z)  ' 

V(y,z)  =  -V(-y,z)  >  (2) 

W(y,z)  =  -W(y,-z)  , 
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and 

fjU(y.O)  -fjV(y,0)  -|^H(0,z)  -  0 

The  displacement  field  of  equations  (1)  should  satisfy  the  equations  of 
equilibrium  and  the  stress-free  conditions  on  the  edges,  y  =  ±b,  and  top  and 
bottom  surfaces,  z  =  ±2h. 

Because  exact  solutions  are  not  available  for  the  complex  three-dimensional 
problem,  approximate  methods  were  used  in  references  1  through  3.  The  present 
study  was  based  on  a  finite-element  analysis,  like  that  in  reference  3,  except 
eight -noded  quadrilateral  isoparametric  elements  were  used  rather  than  three- 
noded  triangular  elements. 

Finite-Element  Idealizations 

Because  of  the  symmetries  in  the  problem,  only  one  quadrant  of  an 
x  =  constant  plane  was  idealized.  The  quadrant  (0  <  y  <  b  and  0  <  z  <  2h) 
is  shown  as  the  shaded  region  in  figure  1(b).  The  displacement  functions  U 
and  V  were  prescribed  as  zero  on  the  y  =  0  line,  and  the  displacement 
function  W  was  prescribed  as  zero  on  the  z  =  0  line. 

The  shaded  region  in  figure  1(b)  was  idealized  by  eight-noded  isoparametric 
elements  in  the  rectangular-mesh  models  shown  in  figure  2  and  the  polar-mesh 
model  in  figure  3.  The  rectangular-mesh  model  had  a  rectangular  arrangement 
throughout.  In  contrast,  the  polar  model  had  a  polar  mesh  near  the  free  edge 
and  rectangular  arrangement  elsewhere.  The  rectangular  mesh  was  used  to  com¬ 
pute  the  stress  distributions,  and  the  polar  mesh  was  used  to  investigate  the 
stress  singularities. 

To  study  the  convergence  of  the  stresses  near  the  free  edge,  three  rec¬ 
tangular  meshes  were  used.  The  medium  mesh  in  figure  2(b)  was  obtained  by 
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subdividing  each  element  of  the  coarse  mesh  (fig*  2(a))  into  four  elements. 
Similarly,  the  fine  mesh  in  figure  2(c)  was  obtained  by  subdividing  each  ele¬ 
ment  of  the  medium  mesh  into  four  elements.  The  coarse  mesh  had  135  nodes  and 
36  elements,  the  medium  mesh  had  485  nodes  and  144  elements,  and  the  fine  mesh 
had  1833  nodes  and  576  elements. 

The  polar  mesh  had  the  polar  arrangement  in  the  region  (19h  <  y  <  20h, 

0  <  z  <  2h)  near  the  laminate  free  edge.  As  in  the  rectangular-mesh  case, 
three  meshes — coarse,  medium,  and  fine — were  used.  Only  the  polar  fine  mesh 
is  shown  in  figure  3.  The  three  polar  meshes  had  the  same  number  of  nodes  and 
elements  as  their  corresponding  rectangular  counterparts.  Near  the  free  edge, 
however,  the  polar  fine  mesh  was  an  order  of  magnitude  finer  than  the  rectangu¬ 
lar  fine  mesh. 

Analysis  of  Singularities 

The  stress  distribution  along  a  radial  line  from  a  singularity  can  be 
expressed  as  (refs.  4  and  5) 

a  -  Alr-°  +  o(r-“+1)  (3) 

where  r  is  the  distance  from  the  singularity,  and  a  are  constants, 

and  0  (r'“+1)  represents  terras  of  the  order  r  0+1  and  higher.  For  small 
distances  r,  the  singular  term  dominates  and  equation  (3)  can  be  approxi¬ 
mated  by 

.  -a 
o  -  A^r 

or 

log  O  =  log  A1  - 


a  log  r 


(4) 
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Hence,  a  log  G  vs.  log  r  plot  would  be  a  line  with  a  slope  of  -a  and  a 
0— intercept  of  A^.  The  d  and  are  the  power  and  the  strength  of  the 

singularity,  respectively. 

Equation  (4)  provided  the  basis  for  a  precedure  to  investigate  the 
existence  of  stress  singularities  at  the  laminate  edge.  This  equation  was 
fitted  to  the  computed  stresses.  If  the  equation  fitted  these  stresses  well 
and  if  the  slope  of  the  fitted  line  was  negative,  a  stress  singularity  with 
power  ot  was  indicated.  This  log-linear  procedure  was  verified  in  refer¬ 
ence  4  by  applying  it  to  several  well-known  singular-stress  problems  in  two- 
dimensional  elasticity. 

The  log— linear  procedure  works  best  where  the  singular  term  strongly 
dominates  the  other  terms.  Therefore,  the  procedure  should  be  applied  only  in 
the  regions  very  close  to  the  singularity. 

RESULTS  AND  DISCUSSION 

First,  distributions  are  presented  for  interlaminar  normal  stress  O z 

and  the  shear  stress  0  for  the  various  [0/ (0-90) ]  -  0  <  0  <  90,  laminates. 

xz  13 

Then,  the  results  are  investigated  for  the  existence  of  stress  singularities 
at  the  free  edge.  For  convenience,  the  intersection  of  the  interface  and  the 
free  edge  (z  =  h;  y  =  b)  will  be  called  the  interface  corner.  Similarly,  the 
intersection  of  the  midplane  and  the  free  edge  (z  =  0;  y  =  b)  will  be  called 
the  midplane  corner . 

The  uniform  axial  strain,  eq,  was  arbitrarily  set  equal  to  0.001  through¬ 
out  this  study. 

Stress  Distributions 

The  stresses  with  the  steepest  gradients  near  the  free  edge  were  the 

normal  stress  a  and  the  shear  stress  O  .  The  distributions  for  these 
z 


8 


stresses,  through  the  thickness  and  along  the  interface,  are  presented  and  com- 

pared  in  this  section. 

Interlaminar  normal  stress,  G  . - 

_  z 

[0/90]  laminate:  Figure  4(a)  shows  the  G  distribution  through  the 

S  2 

thickness  along  the  free  edge  (y  =  b)  for  the  three  rectangular  meshes  (coarse, 

medium,  and  fine).  In  this  figure  and  all  subsequent  figures,  tire  fine  mesh 

results  are  represented  by  a  curve  through  the  data;  only  the  value  at  the 

interface  is  shown  as  a  discrete  value  (diamond  symbol) .  The  coarse  and  medium 

mesh  results  are  shown  by  circular  and  square  symbols,  respectively.  The  solid 

symbols  indicate  the  stresses  in  the  top  ply.  As  shown  on  the  figure,  the 

values  of  G  for  the  three  meshes  agree  closely  except  near  the  interface 
z 

(z  =  h) .  At  the  interface  the  three  meshes  produced  noticeably  different  values. 

In  addition,  G  is  discontinuous  across  the  interface;  this  discontinuity  is 
z 

slightl}?  larger  with  progressive  mesh  refinement.  Hence,  the  stress  state  shown 
in  figure  4(a)  suggests  that  a  singularity  exists  at  the  interface  corner.  In 
contrast,  at  the  midplane  (z  =  o)  the  G^  values  agree  for  the  three  meshes 
and  show  no  evidence  of  a  singularity  at  the  midplane  corner. 

Figure  4(b)  shows  the  average  G^  from  both  plies  plotted  against  a  nor¬ 
malized  distance  from  the  edge.  The  G^  results  from  the  three  meshes  are  in 
excellent  agreement  for  (b  -  y)/h  >  0.08.  However,  at  the  free  edge,  y  =  b, 
the  computed  G^  values  are  noticeably  different  for  the  three  meshes,  with 
the  fine  mesh  producing  the  largest  value.  Therefore,  results  in  both 
figures  4(a)  and  4(b)  suggest  that  a  singularity  exists  at  the  interface  corner. 

[ ±4 5 ]  laminate:  Figure  5(a)  shows  the  G  distribution  through  the 
s  z 

thickness  at  the  free  edge  (y  =  b) .  As  in  the  [0/90]  laminate,  the  values  of 

a  for  the  three  meshes  agree  closely  everywhere  except  near  the  interface 
z 

(z  =  h),  where  a  steep  gradient  exists. 


Figure  5(b)  shows  the  distribution  of 
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O  along  the  interface.  As  expected,  the  results  from  the  three  meshes  agree 
z 

except  very  near  the  free  edge.  Further,  the  computed  value  of  Oz  at  the 

free  edge  was  larger  for  progressive  mesh  refinement.  These  results  for  the 

[+45]  laminate  also  suggest  a  singularity  at  the  interface  corner. 

Other  laminates:  Figures  6  and  7  show  the  O  distributions  through  the 

thickness  and  along  the  interface  for  [15/“75]g,  [75/-15]g,  [30/-60]g,  and 

[ 60 /  —30]  ,  respectively.  In  these  figures  only  the  fine  mesh  results  are 

shown.  The  distributions  of  O ^  near  the  free  edge  have  steep  gradients  for 

all  the  four  laminates,  indicating  a  possible  singularity. 

Comparisons:  The  results  in  figures  4  through  7  show  that,  except  for 

[0/90]  and  [90/0]  laminates,  the  interlaminar  normal  stress,  O ^  is  compres- 
s  s 

sive  at  the  interface  corner.  In  contrast,  at  the  midplane  corner  only  lami¬ 
nates  [75/ —15]  and  [90/0]  develop  compressive  O  . 
s  s 

Interlaminar  shear  stress,  CJXZ*~ 

[±45]  laminate:  Figure  8(a)  presents  the  through-the-thickness  distribu- 
s 

tions  of  a  at  the  free  edge  obtained  using  the  three  rectangular  meshes  for 
xz 

the  [±45]  laminate.  Figure  8(b)  shows  the  distribution  along  the  inter- 

s 

face.  For  all  the  meshes  closely  agrees  everywhere  except  near  the 

interface  in  figure  8(a)  and  free  edge  in  figure  8(b).  The  gradients  are  steep 
in  this  region,  again  indicating  a  possible  singularity  for  the  [±45] g  laminate. 

Other  laminates:  For  the  other  laminates  considered  the  distributions 

were  similar  to  the  [±45] s  results  and,  therefore,  are  not  shown.  Also,  the 
average  0  values  at  the  interface  for  a  laminate  and  its  corresponding  com¬ 
plementary  laminate  (for  example,  [15/-75]g  and  [75/-15]g')  differed  only  in  the 
third  significant  digit. 
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Comparisons:  For  all  laminates  considered  the  0  distributions  had 

the  same  sign  along  the  interface.  The  [ 15/— 7 5]  and  [ 7 5/— 15]  developed  the 

s  s 

largest  0  very  near  the  free  edge  compared  to  other  laminates.  Near  the 

free  edge,  the  0  stresses  are  larger  (by  at  least  a  factor  of  3)  than  the 
xz 

corresponding  stresses  in  each  laminate,  except  for  the  [0/90] g  and  [90/0]g 

laminates  for  which  O  is  identically  zero. 

xz 

For  all  laminates  considered,  O  is  identically  zero  at  the  midplane 

X  Z 

because  of  symmetry.  Therefore,  analysis  of  singularities  is  confined  to  the 
neighborhood  of  the  interface  corner. 

Singularities 

Both  the  0  and  0  distributions  suggest  that  singularities  exist  at 

2  XZ 

the  interface  corner.  To  identify  these  singularities,  both  0  and  0 

stress  data  may  be  used  in  the  log-linear  procedure  outlined  in  the  analysis 

section.  However,  for  the  present  edge-stress  case  (generalized  plane  strain), 

if  one  stress  is  singular  at  a  point,  then  all  stresses  are  and  they  all  have 

the  same  singularity  power,  CL.  This  conclusion  follows  from  the  observation 

that  these  stresses  have  similar  functional  forms,  derived  from  the  same  order 

derivatives  of  the  same  stress  functions,  as  shown  in  reference  6.  Therefore, 

either  0  or  a  may  be  used  in  the  log-linear  procedure  to  determine  the 
2  xz 

power  of  the  singularity. 

To  decide  which  stress  was  better  suited  for  the  log-linear  procedure, 

known  singularities  were  examined.  The  case  of  a  center-cracked  plate  showed 

that  the  log-linear  procedure  gave  best  results  when  the  predominant  stress  was 

used.  For  the  present  study,  the  0  stress  was  predominant.  As  previously 

mentioned,  the  O  values  near  the  edge  were  at  least  three  times  as  large  as 
xz 

the  0  values,  except  for  the  [0/90]  and  [90/0]  laminates.  Further,  as 
z  s  s 


11 


pointed  out  earlier,  the  procedure  works  best  when  the  singular  term  in 
equation  (3)  strongly  dominates  the  other  terms.  The  0^z  distributions  had 
the  same  sign  all  along  the  interface;  in  contrast,  the  0 distributions 
changed  sign  near  the  edge.  This  suggests  that  the  singular  term  for  0^ 
is  dominant  over  a  wider  range  of  data  than  for  .  For  these  reasons 
was  used  in  the  log-linear  procedure. 

The  results  for  the  rectangular  and  polar  meshes  were  indistinguishable 

for  distances  greater  than  0.04h  from  the  free  edge.  Very  near  the  free 

edge,  however,  the  polar  mesh  yielded  a  greater  concentration  of  data  and, 

therefore,  was  used  to  analyze  stress  singularities. 

[±45]g  laminate.-  Figure  9  shows  the  log-log  plots  of  the  absolute  value 

of  O  ,  along  the  interface  for  a  [ ±45]  laminate,  plotted  against  r.  The 
xz  s 

stresses  obtained  from  the  three  polar  meshes  are  plotted  in  this  figure.  As 
expected,  all  three  meshes  show  a  linear  portion  with  a  negative  slope,  indi¬ 
cating  a  singularity.  Before  calculating  the  power  of  the  singularity,  the 
linear  region  for  each  set  of  data  was  estimated  visually.  (The  two  points 
closest  to  the  edge  were  neglected  because  they  represent  nodes  for  the  ele¬ 
ment  at  the  singularity,  and  this  element  cannot  be  expected  to  yield  accurate 
stresses.)  Stresses  from  the  coarse,  medium,  and  fine  meshes  appear  to  be 
linear  in  the  ranges  of  0.0125  <  r  <  0.07,  0.0035  <  r  <  0.015,  and 

0.0015  <  r  <  0.01,  respectively.  In  these  ranges  a  linear  least-square¬ 
fitting  procedure  was  used  to  determine  the  power  (ct)  and  the  strength  of 

the  singularity  in  the  equation 

log  0  “  log  A  -  a  log  r 

XZ  J- 


The  following  values  were  found  for  the  three  meshes. 
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POLAR  MESH 

COARSE 

MEDIUM 

FINE 

Mesh  index,  i 

1 

2 

3 

Power,  a 

0.256 

0.182 

0.170  • 

Percent  change  in  0, 

a.  -  a 

— - —  (100) 

a. 

i 

- 

41 

7 

Strength  A^,  MPa 

5.702 

7.509 

7.950 

As  shown  in  the  table,  by  refining  the  coarse  mesh  to  generate  the  medium 
mesh,  the  a  changed  by  41  percent.  For  the  fine  mesh,  the  additional  refine¬ 
ment  produced  only  a  7— percent  change,  which  is  only  about  one-sixth  the 
42-percent  change.  Further  mesh  refinement  (dividing  each  fine  mesh  element 
into  four  elements)  would  probably  cause  an  additional  1— percent  change 
(one-sixth  of  7  percent).  This  1-percent  change  is  neglibibly  small.  Similar 
arguments  can  be  made  for  the  strength  of  the  singularity,  A±.  Therefore, 
the  present  fine  mesh  results  were  treated  as  converged  values. 

To  check  the  accuracy  of  the  present  procedure,  the  [±45]  singularity 
was  computed  using  the  stresses  along  nine  different  radial  lines  of  the 

polar  mesh. '  The  calculated  power  of  the  singularity  a  should  be  the  same 
for  all  such  lines.  Figure  10  shows  a  and  the  corresponding  values 

plotted  against  (}),  the  angular  position  of  each  radial  line.  These  values  of 
a  vary  slightly,  from  0.195  at  <f>  =  0  to  an  interface  value  of  0.170  at 
(p  =  tt/2.  At  the  interface,  the  0 


stresses  were  largest,  as  shown  by  the 
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largest  A1  value  at  <f>  =  tt/2  in  figure  10.  As  mentioned  earlier,  the  pres¬ 
ent  procedure  yields  best  results  when  the  predominant  stresses  are  used. 

Because  the  O  stress  was  largest  at  the  interface  (<p  =  tt/2),  the  a  value 
xz 

of  0.170  is  believed  to  be  accurate.  Also,  the  absolute  value  of  the  correla¬ 
tion  coefficient  in  the  least— square  procedure  was  closest  to  unity  (0.9996) 
for  <)>  =  tt/2. 

Other  laminates.-  The  O  stress  was  also  used  to  identify  the  singu- 

- - —  XZ 

larities  for  the  [15/-75]  ,  [75/-15]  ,  [30/-60]  and  [60/-30]  laminates. 

s  S 

Figure  11  shows  the  log-log  plots  for  these  four  laminates.  For  comparison 

the  previous  results  for  the  [i45]g  laminate  are  also  included  in  this  figure. 

As  pointed  out  earlier,  for  a  laminate  and  its  complementary  laminate 

the  a  results  differed  only  in  the  third  significant  digit.  Therefore, 
xz 

the  power  of  the  singularity  for  a  laminate  and  its  complement  is  apparently 
identical.  Also  note  that  log-linear  plots  in  figure  11  are  parallel  to  each 
other.  Therefore,  the  power  of  the  singularities  for  [l5/-75]g,  [75/— 15] 

[30/-60] g ,  and  [60/-30] g  are  identical  to  that  for  the  [±45]g  case,  that  is, 
0.170.  However,  the  three  parallel  lines  in  figure  11  are  separated,  indicating 
different  singularity  strengths,  A^ 

Figure  12  shows  a  and  A^  for  various  [0/(0-9O)]g  laminates  as  functions 
of  the  outer  ply  angle,  0.  To  investigate  a  and  A_L  in  the  region 
0  <  0  <  15  (and  75  <  0  <  90),  the  laminates  [2/— 88] g  and  [7.5/-82.5]g  were 
also  considered.  The  a  values  for  these  laminates  are  nearly  equal  to 
0.170,  but  the  A^  values  are  much  smaller  than  for  the  other  laminates.  This 
was  expected  because  the  o xz  stress  is  identically  zero  for  [0/90] g  and 
[90/0]  laminates;  thus,  for  0  nearly  equal  to  zero  A1  should  also  be 
very  small.  Dotted  lines  in  figure  12  represent  extrapolations  to  zero. 
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Figure  12  shows  that  both  0 1  and  A^  are  symmetric  about  0  45°.  Also, 

the  strength,  A^,  is  maximum  for  the  [15/-75]s  and  E 7 5/  1 5 ] ^  laminates. 

[0/90]  and  [90/0]  laminates.-  For  the  [0/90]  laminate,  a  was  used 
_ _s _ _ s _ _ _____  s  z 

in  the  log-linear  procedure  because  the  G  was  identically  zero.  As 

xz 

expected,  the  log-log  plot  in  figure  13  shows  a  linear  region  with  a  negative 
slope.  The  power  of  the  singularity  was  obtained  from  the  least-square  proce¬ 
dure  as  0.202.  This  power  is  larger  than  the  0.17  calculated  from  other 
laminates  but  is  only  slightly  beyond  the  range  shown  in  figure  10.  The  com¬ 
puted  power  of  0.202  probably  occurred  because  the  G^  singularity  for  the 

[0/90]  laminate  was  not  as  dominant  as  the  O  singularity  analyzed  for  the 
s  xz 

other  laminates.  However,  if  the  log-linear  procedure  were  applied  to 
stresses  for  smaller  r  values,  the  power  may  also  approach  0.17.  Such  cal¬ 
culations  for  smaller  r  values  would  require  further  mesh  refinement,  which 
is  beyond  the  scope  of  this  paper. 

For  the  [90/0]  laminate  neither  a  nor  G  could  be  used  in  the 
s  z  xz 

loc— linear  procedure.  The  G  was  identically  zero  and  the  G  changed 
°  xz  z 

sign  very  near  the  free  edge  (see  fig-  4(b)).  Consequently,  the  singularity 
for  the  [90/0]  laminate  was  not  computed.  However,  as  shown  previously  for 
other  laminates,  a  laminate  and  its  complementary  laminate  have  the  same  a.  . 
Therefore,  by  induction,  the  singularities  for  [0/90]  and  [90/0] g  should  also 
be  identical  and  a  should  be  0.17. 

In  general,  these  results  suggest  that  the  power  of  the  singularities  for 

[0/ (0—90) ]  ,  0  <  0  <  90,  laminates  with  material  properties  used  in  this  paper 

s 

is  about  0.17.  This  power  should  be  interpreted  only  as  an  approximate  value 
because  of  the  numerical  accuracy  of  the  log-linear  procedure. 
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CONCLUDING  REMARKS 

The  stresses  near  a  straight,  free  edge  in  composite  laminates  subjected 

to  uniform  axial  strain  were  studied  by  a  quasi-three-dimensional  finite- 

element  analysis.  Seven  different  four-ply  laminates  were  analyzed:  [0/90]  , 

s 

[15/-75]  ,  [30/-60]  ,  [ ±4 5]  ,  [60/-30]  ,  [75/-15]  ,  and  [90/0]  —belonging  to 
s  s  s  s  s  s 

the  family  [0/ (0—90) ]  ,  0  <  6  <  90.  The  laminates  were  idealized  by  eight- 

noded  isoparametric  elements.  Systematic  convergence  studies  were  made  with 
the  aid  of  three  meshes  to  examine  the  interlaminar  stresses  near  the  free 
edge. 

Near  the  intersection  of  the  interface  and  the  free  edge,  the  interlaminar 

normal  stress  o  was  tensile  for  [0/90]  and  [90/0]  laminates  and  was  com- 
z  s  s 

pressive  for  all  other  laminates  considered.  The  interlaminar  shear  stress 

O  was  much  larger  than  the  normal  stress  O  for  all  laminates  except  for 
xz  z 

[0/90]  and  [90/0]  laminates  (for  which  O  is  identically  zero).  Among 
s  s  xz 

the  laminates  considered,  laminates  [15/  — 75]  and  [ 7 5/  15]  developed  the 

s  fa 

highest  O 

XZ 

Convergence  studies  indicated  that  stress  singularities  probably  exist  at 
the  intersection  of  the  interface  and  the  free  edge,  and  the  singularities  do 
not  exist  at  the  midplane.  A  log-linear  curve  fitting  procedure  appears  to 
confirm  the  existence  of  a  singular  point  at  the  intersection  of  the  interface 
and  free  edge.  The  shear  stress  data  close  to  the  singular  point  were  used  in 
the  log-linear  procedure  to  evaluate  the  power  of  the  singularities. 

The  present  results  suggest  that  the  power  of  the  singularity  is  the  same  for 
all  orthotropic  laminates  in  the  family  represented  by  [0/(8-9O)]g,  0  <  0  <  90 

for  a  given  set  of  material  properties.  For  these  laminates,  with  material  pro¬ 
perties  of  the  graphite/ epoxy  considered,  the  power  of  the  singularity  was  about 


0.17. 
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Figure  1.-  Laminate  configuration,  loading,  and  stresses. 
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